Thermodynamics and Statistical Mechanics of frozen systems in inherent states 
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We discuss a Statistical Mechanics approach in the manner of Edwards to the "inherent states" 
(defined as the stable configurations in the potential energy landscape) of glassy systems and granular 
materials. We show that at stationarity the inherent states are distributed according a generalized 
Gibbs measure obtained assuming the validity of the principle of maximum entropy, under suitable 
constraints. In particular we consider three lattice models (a diluted Spin Glass, a monodisperse 
hard-sphere system under gravity and a hard-sphere binary mixture under gravity) undergoing a 
schematic "tap dynamics", showing via Monte Carlo calculations that the time average of macro- 
scopic quantities over the tap dynamics and over such a generalized distribution coincide. We also 
discuss about the general validity of this approach to non thermal systems. 



I. INTRODUCTION 

There are many complex systems where thermal fluc- 
tuations are small enough that the temperature of the 
external bath, Tb at h, can be considered zero. Examples 
are supercooled liquids quenched at zero temperature in 
metastable states (blocked configurations) , called "inher- 
ent structures" , which correspond to the local minima of 
the potential energy in the particles configuration space 
Jl], [|, H . Granular materials || at rest are another im- 
portant example of system frozen Q in mechanically sta- 
ble microstates (blocked configurations), which by anal- 
ogy with the glass terminology can also be called inherent 
states. 

The issue we consider here, which recently raised con- 
siderable interest, is to investigate the possibility to de- 
scribe these systems by using concepts from Statistical 
Mechanics, as Edwards suggested for granular media 
more than ten years ago. His assumption was that, by 
gently shaking the system under the constrain of fixed 
volume V, the distribution over the mechanically sta- 
ble (blocked) states would be uniform. This leads to 
the definition of a configurational entropy, S = lnfi, 
where f2 is the number of mechanically stable states cor- 
responding to the fixed volume V and energy E, and to 
the concept of compactivity, X~ x = d ]^J l ■ In a similar 
way one can also define a configurational temperature, 



a inn 



conf — yconj — q e 

Also in glasses, following for example the inherent 
structure approach [jl], |] |[ 0], one can define a config- 
urational entropy associated to the number of inherent 
structures corresponding to a fixed energy, E, and conse- 
quently the configurational temperature. When the sys- 
tem is frozen at zero temperature in one of his inherent 
states it does not evolve anymore. However, one can 
explore the inherent structures space essentially in two 
ways. One way is by quenching the system over and over 
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from an equilibrium temperature, T, to zero tempera- 
ture (I], ||, |j. Using this procedure, Sciortino et al. || 
found that in a supercooled glass forming liquid, studied 
by molecular dynamics simulations, the configurational 
temperature numerically coincides with the equilibrium 
temperature T, provided that T is low enough. 

Another way to visit the inherent structures is by let- 
ting the system aging in contact with an almost zero 
bath temperature, T^ath- During the aging process an 
effective temperature, Td yn , can be defined via the off- 
equilibrium extension of the Fluctuation-Dissipation ra- 
tio B. It happens that in mean field models this 
effective temperature coincides with the above configu- 
rational temperature. The possibility to introduce an ef- 
fective temperature for granular media via the extension 
of the Fluctuation-Dissipation relation, was suggested in 
Ref. |§. 

The connection between Edwards approach for granu- 
lar media and the results in glass theory has been pointed 
out in 0, [D], [L2| and In particular in Ref. Q it 

was shown that, for a class of finite dimensional systems, 
in the limit Tbath — ► 0, Tdyn coincides in fact with the 
configurational temperature, predicted by the Edwards' 
hypothesis. 

In Ref. jl3] |l4j the inherent states are visited in an 
other way by using a tap dynamics (i.e., a procedure 
similar to that used in the compaction of real granular 
materials), where each tap consists in raising the bath 
temperature to a value Tp and, after a lapse of time To, 
quenching it back to zero. By cyclically repeating the 
process the system explores the space of the inherent 
states Qg^g^ggpOl. Once the sta- 
tionary state is reached one can define a temperature, 
Tfd, via the equilibrium Fluctuation-Dissipation relation. 
One can then see that, if Edwards' assumption applies, 
Tfd coincides with the configurational temperature. This 
has been verified in fact for different finite dimensional 
models g 0, |l5f - It was also shown numerically that 
for low enough Tp one has that Tfd = T con f ~ Tp, con- 
firming on lattice models for granular media the result 
of Ref. J3|. In fact when the duration of each single 
tap is infinite (tq — > oo), the tap coincides with the way 
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to explore the inherent states implemented in molecular 
dynamics simulations for Lennard Jones mixtures (l], [|. 
However the method used in || only allows the calcu- 
lation of T con f when the configurational temperatures is 
low, i.e., where all the different temperatures almost coin- 
cide. Many other studies confirming Edwards' approach 
have also been presented |l6|, [l8|, ^lj . 

In this paper we g ive a comprehensive view of the re- 
sults obtained in fl3[ 0, by considering other models 
and giving more details. In particular we study here three 
schematic lattice models for glassy systems and granular 
media, i.e., a diluted Spin Glass, a monodisperse hard- 
sphere system under gravity and a hard-sphere binary 
mixture under gravity. In particular, in the diluted Spin 
Glass and in the monodisperse hard-sphere system un- 
der gravity, the asymptotic states reached by the system 
are found to be described only by the configurational 
temperature. Whereas in the hard-sphere binary mix- 
ture under gravity the asymptotic states are found to be 
described by two thermodynamic parameters [p2[ , coin- 
ciding with the two configurational temperatures which 
characterize the distribution among the inherent states 
when the principle of maximum entropy is satisfied un- 
der the constraint that the energies of the two species are 
independently fixed. In Ref. |H| a description of the seg- 
regation observed in the binary system in term of these 
two temperatures is also given. 



In Sect.s [I A and II B , the frustrated lattice gas model 
and the results of its study with the tap dynamics are 
respectively presented. In Sect. |EI D| , the same results 
are obtained at higher density where the system at small 
temperature reaches a quasi-stationary state in which 
one ti me q uantities decay as the logarithm of time. In 
Sect. [IC Edwards' hypothesis is formulated using the 
principle of maximum entropy. The results obtained in 
the monodisperse hard-sphere system under gravity are 
shown in Sect. [II. In Sect. |[y|the Statistical Mechanics 
approach is extended to the hard-sphere binary mixture 
under gravity, where two thermodynamic parameters are 
necessary to describe the asymptotic states reached by 
the system. Finally, in the Conclusions we draw a picture 
of the Statistical Mechanics approach to systems found 
in inherent states, as it emerges from our extensive in- 
vestigation. 



where the sum ^ s over nearest neighbor sites, Si — 

±1 are Ising spins, rii = 0, 1 are occupation variables, 
fi is the particle chemical potential and are quenched 
and random variables, equal to ±1 with equal probability. 
This model reproduces the Ising spin glass in the limit 
p — ► oo (i.e., when all sites are occupied, rii = 1). 

In the other limit, J — > oo, the model describes a 
frustrated lattice gas with properties recalling those of a 
"frustrated" liquid. In fact the first term of Hamiltonian 
(|l|) implies that two nearest neighbor sites can be freely 
occupied only if their spin variables satisfy the interac- 
tion, that is if CijSiSj = 1, otherwise they feel a strong 
repulsion. To make the connection with a liquid, we note 
that the internal degree of freedom, Si, may represent 
for example the internal orientation of a non spherical 
particle. Two particles can be nearest neighbors only 
if the relative orientation is appropriate, otherwise they 
have to move apart. Since in a frustrated loop the spins 
cannot satisfy all interactions, in this model particle con- 
figurations in which a frustrated loop is fully occupied are 
not allowed. The frustrated loops in the model are the 
same of the spin glass model and correspond in the liquid 
to those loops which, due to geometrical hindrance, can- 
not be fully occupied by the particles. In 31? |||, Hg], the 
model has a maximum density /o max — 0.68, and a transi- 
tion at p c ~ 0.62 where the non-linear spin susceptibility 
diverges. 

In the present paper, the 3D cubic frustrated lattice 
gas model with J finite is considered. The value of parti- 
cle density, p = rii/L 3 (L is the lattice linear size), is 
fixed, and a Monte Carlo tap dynamics, which allows the 
system to explore its inherent states, is applied. During 
the dynamics, the system cyclically evolves for a time To 
(the tap duration |p0| ) at a finite value of the bath tem- 
perature, Tr (the tap amplitude), and afterwards it is 
suddenly frozen at zero temperature in one of its inherent 
states (at zero temperature the system does not evolve 
anymore if the energy cannot be decreased by one single 
particle movement). After each tap, when the system is 
at rest, we record the quantities of interest. The time, 
t, considered is therefore discrete and coincides with the 
number of taps. 



B. The results obtained under the tap dynamics 



II. THE FRUSTRATED LATTICE GAS MODEL 



The model 



Recently a lattice model has been introduced to de- 
scribe glass formers [p5, |26| a nd, in presence of gravity, 
granular materials Jl9|, 27|, B8j. The Hamiltonian of the 
model is: 



— H = J^^(eijSiSj — l)nirij 



(1) 



(ij) 



We first consider the case p — 0.65 j3lj . At this value of 
the density the system under the tap dynamics reaches a 
stationary state for each value of Tr (and to) considered. 
In Fig.s [l]and[| |32| , the self-scattering two time function, 
F q (t,t w ) = Y^eFF'W-^^/pL 3 , and the energy, E(t), 
of the inherent states, obtained for Tr = 0.3 J and tq — 
1 MCS, are shown. The curves, F q (t,t w ), for different 
t Wl collapse onto a single master function, when they 
are plotted as function of t — t Wl and the energy, E(t), 
reaches its time independent asymptotic value, showing 
that the system has reached a stationary state (our data 
are averaged up to 32 noise realizations; L = 8 and q = 
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FIG. 1: The self-scattering two time function, F q (t,t w ) ~ 
Y^.e^^-^^/pL 3 , with <? = tt/4, as function of t - t w 
(for t w = 10 4 , 2-10 4 , 5-10 4 , 8-10 4 , 10 s ) in the frustrated lattice 
gas model for density p = 0.65, during the tap dynamics, with 
tap amplitude Tr = 0.3 J and tap duration to — 1 MCS. The 
function F q (t,t w ) only depends on t — t m , showing that the 
system has reached stationarity. 



FIG. 2: The energy, E(t), of the inherent states as function of 
the tap number t, in the frustrated lattice gas model during 
the tap dynamics with tap amplitude Tr = 0.3 J and tap 
duration To = 1 MCS. The lower curve, corresponding to a 
density of particles p — 0.65, exponentially saturates to its 
asymptotic value, whereas the upper curve, corresponding to 
p = 0.75, shows a logarithmic relaxation at long times. 



tt/4). 

During the tap dynamics, in the stationary state, the 
time average of the energy, E, and its flu ctuat ions, AE 2 , 
are calculated. In Fig.s || and ||, E and AE 2 are shown 
as function the tap amplitude, Tr, (for several values of 
the tap duration, To). Apparently, Tr is not the right 
thermodynamic parameter, since sequences of taps, with 
same T r and different To, give different values of E and 
AE 2 . On the other hand, if the stationary states are 
indeed characterized by a single thermodynamic param- 
eter, f3fd, the curves corresponding to different tap se- 
quences (i.e. different Tr and To) should collapse onto a 
single master function, when AE 2 is parametrically plot- 
ted as function of E. This data collapse is in fact found 
and shown in Fig. || This is a prediction which could be 
easily checked in real granular materials (where one could 
consider the density which is easier to measure than the 
energy). 

The thermodynamic parameter, [3fd, is defined apart 
from an integration constant, (3q, through the usual equi- 
librium Fluctuation-Dissipation relation: 



C. The Edwards' averages 



dE 



AE 2 . 



(2) 



By integrating eq. (0), (3fd — 0o can be expressed as 
function of E or (for a fixed value of To) as function of 
fir- Pfd — A) = 9{Pr)- The constant /?o is determined as 
explained in detail in the Appendix. 



In Sect. II B we have found that the fluctuations of the 



energy in the stationary state depend only on the energy, 
E, and not on the past history. If all macroscopic quan- 
tities depend only on the energy, E, or on its conjugate 
thermodynamic parameter, f3fd, the stationary state can 
be genuinely considered a "thermodynamic state" . In 
this case one can attempt to construct an equilibrium 
statistical mechanics, as originally suggested by Edwards 

More precisely one can try to find from basic gen- 
eral principles what is the probability distribution, P r , 
of finding, in the stationary regime, the system in the 
inherent state r of energy E r (see |H). We assume 
that the distribution is given by the principle of max- 
imum entropy, S — — P r lnP r , under the condi- 
tion that the average energy is fixed: E = ^2 r P r E r . 
Thus, we have to maximize the following functional: 
I[Pr] = -J2 r p r^Pr-Pconf(E-J2 r P r E r ). Here [3 conf 
is a Lagrange multiplier determined by the constraint on 
the energy and takes the name of "inverse configurational 
temperature". This procedure leads to the Gibbs result: 



e -PconfE r 



(3) 



where Z = J2 r e~" aon f Er . Using standard Statistical Me- 
chanics it is easy to show that, in the thermodynamic 
limit, the entropy S and (3 con f are also given by: 



S = ]nQ(E), /3 conf 



d In fl 
~dE~ 



(4) 
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FIG. 3: The time average of the energy, E, recorded in the 
stationary regime, as function of the tap amplitude, Tr (in 
units of J), in the frustrated lattice gas model for p = 0.65. 
The four different curves correspond to different values of the 
tap duration, to = 1, 5, 10, oo MCS (from bottom to top). 
This shows that Tr is not a right thermodynamic parameter, 
since sequences of taps with different to give different values 
for the system observables. 
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FIG. 4: The time average of energy fluctuations, AE 2 , 
recorded in the stationary regime, as function of the tap am- 
plitude, Tr (in units of J), in the frustrated lattice gas model 
for p — 0.65. The four different curves correspond to dif- 
ferent values of the tap duration, to = 1, 5, 10, oo MCS 
(from bottom to top). This shows again that Tr is not a right 
thermodynamic parameter. 



where fl(E) is the number of inherent states correspond- 
ing to energy E. 

If the distribution in the stationary state coincides with 
eq.(||) the time average of the energy, E(/3fd), recorded 
during the taps sequences, must coincide with the en- 
semble average, (E)(/3 con f), over the distribution eq.(||). 
In order to check that we have independently calculated 
the average (E), as function of f3 con f. We have simu- 
lated the model eq.(|l|) imposing that the only accessible 
states are the inherent states, as done in Ref. |0|]. The 
only difference is that in the present paper the Edwards' 
averages are done in the canonical ensemble, whereas 
in Ref. [ pT| are done in the microcanonical ensemble. 
In particular we have constructed a new Hamiltonian, 
H'({Si,ni}) = H({Si,ni}) + *({S<,ni}), by adding a 
term to eq.(Q), 5({Si,rii}), which is zero, if the con- 
figuration is an inherent state, and infinite, otherwise. 
The canonical distribution for this Hamiltonian gives a 
weight, e~' 9co " / ' H , which coincides with the weight in the 
distribution of eq.(||), for each accessible configuration. 
Using the standard Monte Carlo simulations, we have 
calculated (E)(f3 con f). Fig. || outlines a very good agree- 
ment between (E)((3 con f) and E(f3fd) (notice that there 
are no adjustable parameters). 

Note that the maximum energy reached by the system 
under the tap dynamics, E max (To) = E(Tp — > oo,to), 
is less than the maximum energy of the inherent states, 
(E)(T con f — > oo), for every value of r considered. Such 
a prediction, which may have important practical con- 
sequences (e.g., in powders technologies), is consistent 



with some experimental observations on tapped granular 
materials |33j , where the system density was shown to ap- 
proach asymptotically a plateau value apparently higher 
than the minimal possible packing density (obtained, for 
instance, by just pouring grains in the container) even 
for very large tap amplitudes. 

Using eq.([|), we have finally evaluated the configura- 

tional entropy as S(E) — So = f Q (3 con f(E')dE' (where 
the unknown non negative constant, So = S(E — 0), is 
the entropy at T con f = 0). In Fig. |?j, the configurational 
entropy, S — So, is plotted as function of T con f. We have 

also evaluated the integral S'(E) -S' = J® /3 fd (E')dE'. 
In Fig. f?| S' — S'o is plotted as function of Tjd and it is 
compared with the configurational entropy. The agree- 
ment is again very good. 



D. Quasi-stationary case 

We have also studied the frustrated lattice gas model 
for p = 0.75. Differently from the previous case, for 
small enough values of the tap amplitude Tr, the system 
does not reach a stationary state during our observation 
time. In Fig. 0, the energy, E(t), of the inherent states 
obtained for T r = 0.3 J and r = 1 MCS is shown. E(t) 
now changes in time and the system is not in a stationary 
state; however, E(t), at long times decays very slowly 
pij . In this regime the time averages are computed over a 
time interval such that the energy is practically constant 
(in the case of Fig. |[ the time average is performed 
over the time interval (3T0 5 ,3T0 5 + 10 4 )). Performing 
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FIG. 5: The time averages of energy fluctuations, AE 2 , when 
plotted as function of the time average of energy, E, collapse 
onto a single master function, for all the different values of 
tap amplitude and duration, Tr and To, plotted in Fig.[]. This 
shows that the system stationary states are indeed character- 
ized by a single thermodynamic parameter, since the curves 
corresponding to different tap sequences (i.e. different Tr and 
ro) collapse on a "universal" function, when AE 2 is paramet- 
rically plotted as function of E. 
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FIG. 6: The time average, E, calculated in the stationary 
regime of the tap dynamics and the ensemble average over the 
Edwards' distribution eq.(^|), (E) (the black empty circles), 
are plotted respectively as a function of Tfd and T con f (in 
units of J), in the frustrated lattice gas model, at p = 0.65. 
The two independently calculated sets of points show a very 
good agreement, outlining the success of Edwards' approach 
to describe the system macroscopic properties. 



FIG. 7: The configurational entropy, S — So, (the black empty 
circles in figure) as function of T co „/ (in units of J), compared 
with S'{E) - S' = p fd (E')dE' plotted as function of T fd 
(in units of J), in the frustrated lattice gas model for p = 
0.65. The unknown non negative constant So is the entropy 
at T con f = 0. 
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FIG. 8: The time average, E, and the ensemble average over 
the distribution eq.(^), (E) (the black empty circles), plotted 
respectively as a function of Tfd and T con f (in units of J), 
in the frustrated lattice gas model, at p = 0.75. As well as 
at p — 0.65, there is a very good agreement between the two 
independently calculated sets of points. 



the same procedure described in the stationary 
collapse of data is again found (see Fig.||). 

We have again evaluated the configurational entropy, 
S(E) - So = Jz mm (3 C onf(E')dE' (where the unknown 
non negative constant, Sg = S(E min ), is the entropy at 
T con f — and E m i n is the minimum value of energy 
obtained p5|). In Fig. pi, the configurational entropy, 
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FIG. 9: The configurational entropy, S~ So, (the black empty 
circles in figure) as function of T con f (in units of J), compared 
with S'(E) - S' = J E P fd (E')dE' plotted as function of T Jd 
(in units of J), in the frustrated lattice gas model for p — 
0.75. The unknown non negative constant So is the entropy 
at Tconf = 0. 



S — Sq, is plotted as function of T< 



conf • 



Wc have also 



evaluated the integral S'(E) - S' = J Q f3 fd (E')dE'. In 
Fig. ||, 5' — 5q is plotted as function of Tfd and it is com- 
pared with the configurational entropy. The agreement 
is again very good. 

We cannot exclude that the agreement here found even 
for low energy may be due to the fact that the system, 
which is not in a stationary state, is however very near 
stationarity. 



III. A MONODISPERSE HARD-SPHERE 
SYSTEM UNDER GRAVITY 

As a model more appropriate for granular media, we 
have also studied a system of monodisperse hard-sphere 
(with diameter clq = 1) under gravity, where the centers 
of mass of grains are constrained to move on the sites of 
a cubic lattice (see upper inset in Fig. [II]). The Hamil- 
tonian of the system is: 




FIG. 10: The density self-overlap function, Q, and (upper 
inset) the system density on the bottom layer, pb, plotted 
as function of Tfd (in units mgao), compared with the en- 
semble averages over the distribution eq.(p|) (the black empty 
circles), plotted as function of T co „f (in units mgao), in the 
3D monodisperse hard-sphere system under gravity. Also for 
this system there is a very good agreement between the in- 
dependently calculated time averages over the tap dynamics 
and the statistical mechanics ensemble averages a la Edwards. 



to explore its inherent states. We have considered 3 dif- 
ferent values of the tap duration, to = 500, 10, 5 MCS. 
In this case we again obtain that the asymptotic states 
reached by the system can be described by a single ther- 
modynamic parameter, fifd, evaluated by integration of 
eq. (0). We have moreover calculated the system den- 
sity on the bottom layer, pt, and the density self-overlap 
function, Q, and verified that, when plotted as function 
of Pfd, they scale on a single master function (see Fig. 
0). 



H = H h c({n t }) + ffTO i 



(5) 



As described in Sect. [IC, we have calculated the Ed- 
wards' averages as function of (3 C onf- As we can see 
in Fig.|ll|, we obtain a very good agreement between 
(E)(j3 con f) an d E(f3fd). The same agreement is found 
for the other quoted observables, pb and Q (see Fig. [To|). 



IV. A HARD-SPHERE BINARY MIXTURE 
UNDER GRAVITY 



where the height of site i is Zj, g = 1 is the gravity accel- 
eration, m = 1 the grains mass, rii is the usual occupancy 
variable and H.hc{{ni}) is the hard core term preventing 
the overlapping of nearest neighbors grains (the analogy 
with eq.(0) can be appreciated by writing down TLhc'- it 
can be written as Hhc({ni}) = JJ2{ij) n i n ji where the 
limit J — ► oo is taken) . 

We have considered a system of N = 240 particles, 
and performed a tap dynamics which allows the system 



Finally we consider a hard-sphere binary system made 
of two species 1 (small) and 2 (large) with grain diam- 
eters ao and y/2ao, under gravity on a cubic lattice of 
spacing ao = 1. We set the units such that the two kinds 
of grain have masses mi = 1 and = 2, and the grav- 
ity acceleration is g = 1. The hard core potential Hhc 
is such that two large nearest neighbor particles cannot 
overlap. This implies that only couples of small parti- 
cles can be nearest neighbors on the lattice. The system 
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FIG. 11: Main frame The time average, E, and the ensem- 
ble average over the distribution eq.(^), (E) (the black empty 
circles), plotted respectively as a function of Tfd and T con f (in 
units mgao), in the 3D monodisperse hard-sphere system un- 
der gravity described in the text (and schematically depicted 
in the upper inset). Time averages over the tap dynamics 
and Edwards' ensemble averages then coincide. Lower In- 
set The temperature Tfd = /3J^ defined by eq.(Q) as function 
of Tr (in units mgao) for to = 500, 10,5 MCS (from top to 
bottom). The straight line is the function Tfd = Tr. 



overall Hamiltonian is: 



ri = Hhc + rnigHx + m 2 gH 2 



(6) 



where Hi = J2i' 7 -i an( i ^2 = TX' z «> the height of site 



% is Zi and the two sums are over all particles of species 1 
and 2 respectively. In the above units, the gravitational 
energies in a given configuration are thus E± = Hi and 
E-2 = 2H2- 

Grains are confined in a box of linear size L between 
hard walls and periodic boundary conditions in the hori- 
zontal directions. N\ = 120 grains of type 1 and iV2 = 40 
grains of type 2 are initially prepared in a random loose 
stable pack. Under the tap dynamics the system ap- 
proaches a stationary state for each value of the tap pa- 
rameters Tr and tq. In Fig.|l| we plot as function 
of Tr (for several values of to) the asymptotic value of 
the vertical segregation parameter, i.e., the difference of 
the average heights of the small and large grains at sta- 
tionarity, A/i(Tr,To) = h± — h 2 . Here hi and h 2 are 
the average of H\/Ni and H2/N2 over the tap dynam- 
ics in the stationary state. (An interpretation, in terms 
of the approach here presented, of the size segregation 
phenomenon here found and experimentally observed in 
a hard-spheres binary mixture under gravity is given in 
Ref. g§). 

The results given in the main panel of Fig. |l2j appar- 
ently show that Tr is not a right thermodynamic pa- 
rameter, since sequences of taps with different tq give 
different values for the system observables. However, if 




FIG. 12: Main frame The difference of the average heights of 
small and large grains, Ah — h\—h2, measured at stationarity 
in the binary hard spheres mixture under gravity, is plotted 
as a function of tap amplitude, Tr (in units mgao). The 
three sets of points correspond to the shown tap durations, 
To. At high Tr larger grains are found above the smaller, i.e, 
Ah < 0, as in the Brazil nut effect (BNE). Below a Tr(r ) the 
opposite is found (Reverse Brazil nut effect, RBNE). Inset 
The Ah data of the main frame are plotted as a function of 
the corresponding average energy, e. The three sets of data 
do not collapse onto a single master function, showing that a 
single macroscopic observable, such as e, doesn't characterize 
the system status. 
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FIG. 13: The average density of small grains on the box bot- 
tom layer, p 1 , measured at stationarity as a function of the 
height of small particles, hi, in the binary hard spheres mix- 
ture under gravity. Data corresponding to different Tr and to 
approximately scale on a single master function. The empty 
circles are the corresponding values obtained by ensemble av- 
erage with the two temperatures Gibbs measure proposed in 
the text. 
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FIG. 14: Main frame The average density of large grains 
on the box bottom layer, p\, obtained for different Tr and 
ro, scale almost on a single master function when plotted as 
a function of the large grains height, fe. Upper inset The 
average number of contacts between large grains per particle, 
N c , obtained for different Tr and ro, scale on a single master 
function when plotted as a function of the system energy, e. 



the stationary states corresponding to different tap dy- 
namics (i.e., different Tr and To), are indeed character- 
ized by a single thermodynamic parameter, the curves 
of Figjl^ should collapse onto a universal master func- 
tion when Ah(Tr,To) is parametrically plotted as func- 
tion of an other macroscopic observable such as the av- 
erage energy, e(T r ,r ) = (£1 + E 2 )/N (N is the total 
number of particles). This collapse of data is clearly not 
observed here, as it is apparent in the inset of Fig.|l2} 
We show, instead, that two macroscopic quantities may 
be sufficient to characterize uniquely the stationary state 
of the system. These two quantities are, for instance, 
the energy e and the height difference Ah. Of course 
since e = ahi + 2bh 2 (where a = Ni/N and b = N 2 /N) 
and Ah = h\ — h 2 , we can also choose hi and h 2 to 
characterize the stationary state. Namely, we show that 
any macroscopic quantity A, averaged over the tap dy- 
namics in the stationary state, is only dependent on 
hi and h 2 , i.e., A = A(hi,h 2 ). We have checked that 
this is the case for several independent observables, such 
as the number of contacts between large particles, N c , 
the density of small and large particles on the bottom 
layer, and On, and others. In particular, as shown in 
Fig.s and O, we find with good approximation that: 
N c ~ N c (e) = N c (ahi + bh 2 ), p\ ~ p b 2 (h 2 ), p\ ~ p\(hi). 
Therefore we need both hi and h 2 to characterize unam- 
biguously the state of the system; namely all the observ- 
ables assume the same values in a stationary state char- 
acterized by the same values of hi and h 2 , independently 
on the previous history (i.e., in our case independently 
on the particular tapping parameters Tr and To). 

We again find that the stationary state can be gen- 
uinely considered as a thermodynamic state. Therefore 



we can ask what is the probability distribution, P r , of 
finding the system in the inherent state r corresponding 
to an energy Ei r for the small particles and E 2r for the 
large particles. We again assume that the microscopic 
distribution is given by the principle of maximum entropy 
S = — P r In P r , now under the condition that the av- 
erage energy £1 = Y^ r Pr E ir and E 2 = J2 r p rE 2r are in- 
dependently fixed. This can be done by introducing two 
Lagrange multipliers [3\ and (3 2 , which are determined 
by the constraint on E\ and E 2 , and can be respectively 
considered as the "inverse configurational temperature" 
of species 1 and 2. This procedure leads to the Gibbs 
result: 



E\ r — @2 



Pr 



(7) 



where Z = exp \—f3\Ei r — (3 2 E 2r \ and, in the ther- 
modynamic limit, the entropy, S, is given by: 



S = ]nn(Ei,E 2 



and Pi and (3 2 : 



d\nn{E u E 2 ) 
dEi 



P 2 = 



dlnn(gi,g 2 ) 
dE 2 



(8) 



(9) 



Here 17 (£1, £2) is the number of inherent states corre- 
sponding to energy E\ and £2. The hypothesis that the 
ensemble distribution at stationarity is given by eq.(|7]) 
can be tested as follows. We have to check that the time 
average of any quantity, A(hi,h 2 ), as recorded during 
the taps sequences in a stationary state characterized by 
given values hi and h 2 , must coincide with the ensemble 
average, (A) (hi, h 2 ), over the distribution eq.(^). To this 
aim, we have calculated the ensemble averages (N c ), (p 2 ), 
(p\) for different values of (3i and f3 2 ;we have expressed 
parametrically (N c ), (p 2 ), (p\), as function of the average 
of hi and h 2l and compared them with the corresponding 
quantities, A^ c , p\ and p\, averaged over the tap dynam- 
ics. The two sets of data are plotted in Fig.s [l^ and 
|l4| showing a good agreement (notice, there are no ad- 
justable parameters). In order to calculate the ensemble 
averages we simulate the model with Tl from eq.(^) where 
we impose the constraint that the only accessible s tates 
are the inherent states, as already described in Sect. II C . 



V. CONCLUSIONS 

In conclusion, in the context of models for glasses and 
granular materials we have shown that the stationary (or 
quasi-stationary) state reached by the system subject to 
a tap dynamics among its inherent states is genuinely 
a thermodynamic state, which can be well described by 
Edwards' assumption of an uniform measure, i.e., a prob- 
ability distribution obtained assuming the validity of the 
principle of maximum entropy. In particular in the frus- 
trated lattice gas model and in the system of monodis- 
perse hard-spheres under gravity we have found that the 



9 



observables recorded during different tap sequences (dif- 
ferent amplitude and duration of taps) fall on universal 
master curves when plotted as a function of a single ther- 
modynamic parameter. These curves turn out to coin- 
cide with those predicted, within the described Statistical 
Mechanics approach, by the generalized Gibbs distribu- 
tion of eq.(|^). On the other hand, the results obtained 
in a system under gravity made of particles of two dif- 
ferent sizes show that a single thermodynamic parame- 
ter is not enough to describe the macrostates, and two 
configurational temperatures are instead necessary. In 
general, for a more complex system one might expect 
more constraints to be imposed, leading to more than 
two thermodynamic parameters |l5|, |l8|, p3|. In practice, 
the criteria to determine a priori the required parame- 
ters can be not easily accessible. However more recently 
we have extended data regarding the hard-sphere binary 
mixture for very low energy p(| and found that only one 
thcrmodynamical parameter is necessary to describe the 
stationary state. This seems to be a general feature p2). 
If this is the case, a Statistical Mechanics approach with 
only one thermodynamical variable may be feasible for 
low energy. 
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APPENDIX A: DETERMINATION OF THE 
INTEGRATION CONSTANT fa 

Adapting to lattice models the procedure of Sciortino 
et al. ||, we have evaluated (3fd at small values of Tr, 
for To — > oo, and consequently the integration constant, 
ft- 

Given an inherent state, r, of energy, E r , we define the 
basin of attraction, B r , of such state r as the set of states 
in the configurational space, which after the quench at 
T = are frozen in the inherent state r. Therefore the 
probability distribution, P r , of finding the system in the 
inherent state r after quenching the system from an equi- 
librium state at temperature Tr, can be written as 



Pr = 



-■Err' AT 



Zg{Tt) 



(Al) 



where Zq(Ty) is the partition function of the system in 
equilibrium at temperature Tr and ^2 , is the the sum 
over all the states r' belonging to the basin B r of energy 
E rr i . By putting E rr i = E r + A rr /, the distribution (Al) 
can be written as 

e -{E r +g r (T T ))/T T 

Pr = - Z G (T r) ■ ^ 
where e _ff '' (Tr) / Tr = J2 r > e ~ Arr ' • From e q- flU) [t 
follows that the probability of finding the system in any 
inherent state of energy E, P(E,Tr) — ^ r P r (where 
is the sum over all the inherent states r of energy E) 
is given by 



P(E,T T ) 



Q(E)e 



-E/T r -f(T r ,E)/T r 



Zg{Tt) 



(A3) 



where Q(E) is the number of inherent states of energy E 



and e -f(Tv,E)/T T 



'/Tr 



n(E) 



. From eq. (A3) 



\n[P(E,T r ) 



E 



f(T T ,E) 

Tr 



In [Sl(E)]- 



■]n[Z G (T r )]. 
. (A4) 

The probability distribution of finding the system in any 
inherent state of energy E, P{E, Tr), is measured during 
the tap dynamics with amplitude tq — > oo. If f(Tr,E) 
has only a weak dependence on E, then it is possible to 
superimposed the curves, In [P(E, Tr)] + E/Tr, at differ- 
ent Tr which overlap in E by adding a Tr— dependent 
constant. This result is obtained for Tp < 0.525, 
as shown in Fig.na, and suggests that in this i nter val 
f(T r ,E) ~ f(T r y If this is the case, from eq.Q it 
follows: 



P(E,T r ) 



n(E)e 



-E/T T 



Z(T V ) 



where 



Z(T r )=eS^I T -Z G (T r ) 



E 



Sl(E)e 



-E/T r 



(A5) 



(A6) 



The last equality stems from the normalization condition 
on P{E,T T ). 

From eq. flA5|) it follows that at small Tr, /3r = T^ 1 
satisfies eq.(|2|). Therefore at small Tr, f3fd and Pr co- 
incide. The constant /3q is consequently obtained as the 
limit, for Tr — *■ 0, of the function {3t — g(Pr) (see Fig.^6[). 



[1] F.H. Stillinger and T.A. Weber, Phys. Rev. A 25, 978 [2] B. Coluzzi, G. Parisi and P. Verrocchio, Phys. Rev. Lett. 
(1982); Science 225, 983 (1984). F. H. Stillinger, Science, 84, 306 (2000). 

267 1935, (1995). S. Sastry, PG. Debenedetti, F.H. Still- [3] F. Sciortino, W. Kob, P. Tartaglia, Phys. Rev. Lett. 83, 
inger, Nature 393, 554 (1998). 3214 (1999). 



10 



In [P(E,T r )]+T r /E 

90 r 

80 r j 
70 : k * 

60 r . A 

■ 

50 r 

A 

40 : 

30 : f 
20 r / 

10 y° 

n r i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

u 5 10 15 20 25 30 35 40 45, 



FIG. 15: The curves In [P(J3, Tr)] + E/Tr (apart from a 
Tr — dependent constant) as function of the energy E in 
the frustrated lattice gas model for p = 0.65 and Tr = 
0.275, 0.425, 0.475, 0.525 J. 



10 



FIG. 16: The curve (3 T - g(/3 r ) = j3 r - [/8/d(/9r) - A)] as 
function of fi^ 1 (in units of J) in the frustrated lattice gas 
model for p — 0.65. The limit /Jp 1 — > of /3r — g(Pr) gives 
the integration constant /3q- 



[4] W. Kob, F. Sciortino and P. Tartaglia, Europhys. Lett., 
49, 590 (2000). 

F. Sciortino and P. Tartaglia, Phys. Rev. Lett. 86, 107 
(2001). 

[5] H.M. Jaeger, S.R. Nagel and R.P. Behringer, Rev. Mod. 
Phys. 68, 1259 (1996). 

[6] Grains are "frozen" because, due to their large masses 
and dissipation || , the thermal kinetic energy is negligi- 
ble compared to the gravitational energy; thus the exter- 
nal bath temperature, Ttath, can be considered equal to 
zero. 

[7] S.F. Edwards and R.B.S. Oakeshott, Physica A 157, 



1080 (1989). A. Mehta and S.F. Edwards, Physica A 157, 
1091 (1989). S.F. Edwards, in "Disorder in Condensed 
Matter Physics" page. 148, Oxford Science Pubs (1991); 
and in Granular Matter: an interdisciplinary approach, 
(Springer- Verlag, New York, 1994), A. Mehta ed. 

[8] L. F. Cugliandolo, J. Kurchan and L. Peliti, Phys. Rev. 
E 55,3898 (1997). 

[9] R. Monasson, Phys. Rev. Lett. 75, 2847 (1995). Th.M. 
Nieuwenhuizen, Phys. Rev. E 61, 267 (2000). S. Franz 
and M.A. Virasoro, J. Phys. A 33, 891 (2000). A. Crisanti 
and F. Ritort, Journ al of Chemical Phy sics 113, 10615 



[10 
[11 



[12 
[13 
[14 





g(&) 


[15] 


0.7 


-\ 

\ 

\ 


[16] 


0.6 


• 


[17] 


0.5 


m 




0.4 


- m 


[18] 


0.3 


7 • 


[19] 


0.2 


• 








[20] 


0.1 


• 

• 


[21] 





i i 


[22] 



[23 
[24 

[25 

[26 



(2000). J. Kurchan, cond-mat/981Z3J i l\ and in "Jam 



ming and Rheology: Constrained Dynamics on Micro- 
scopic and Macroscopic Scales", A.J. Liu and S.R. Nagel 
Eds., Taylor and Francis, London (2001). 
M. Nicodemi, Phys. Rev. Lett. 82, 3734 (1999). 
A. Barrat, J. Kurchan, V. Loreto and M. Sellitto, Phys. 
Rev. Lett. 85, 5034 (2000); Phys. Rev. E 63,051301 
(2001) 

H. A. Makse and J. Kurchan, Nature 415, 614 (2002). 
A. Coniglio and M. Nicodemi, Physica A 296, 451 (2001). 
A. Coniglio, A. Fierro and M. Nicodemi, Physica A 
302 (1-4): 193 (2001); A. Fierro, M. Nicodemi and A. 
Coniglio, Europhys. Lett. 59, 642 (2002); A. Coniglio, 
A. Fierro and M. Nicodemi, Europ. Phys. Jour. E, in 
press (2002). 

M. Nicodemi, A. Fierro and A. Coniglio, cond- 
maf/0202500 and Europhys. Lett, in press. 
J.J. Brey, A. Prados, B. Sanchez-Rey, Physica A 275, 310 

(2000) ; A. Prados, J.J. Brey, B. Sanchez-Rey, Physica A 
284, 277 (2000). 

D. S. Dean and A. Lefevre, Phys. Rev. Lett. 86, 5639 

(2001) ; J. Phys. A: Math. Gen. 34, L213 (2001); cond- 
roat/0106220. 

J. Berg, S. Franz and M. Sellitto, Eur. Phys. J. B 26, 
349 (2002). 

M. Nicodemi, A. Coniglio, H.J. Herrmann, Phys. Rev. E 
55, 3962 (1997); J. Phys. A 30, L379 (1997). A. Coniglio 
and H.J. Herrmann, Physica A 225, 1 (1996). 
J. Berg, A. Mehta, Europhys. Lett. 56, 784 (2001). 
L. Berthier, L.F. Cugliandolo and J.L. Iguain, Phys. Rev. 
E 63, 051302 (2001). 

the possibility of introducing more than one thermody- 
namic parameter has been also suggested in jl8| and re- 
cently discussed in the context of a Constrained Ising 
Chain in [Q. 

A. Lefevre, cond-mat/0202376. 

A. Mehta and G. C. Barker, Phys. Rev. Lett. 67, 394 
(1991). 

A. Coniglio, A. de Candia, A. Fierro, M. Nicodemi, Jour. 
Phys.: Cond. Mat. 11, A167 (1999). 
A. Coniglio, "Frustration and connectivity in glass form- 
ing systems and granular materials", in Proc. Int. School 
on the Physics of Complex Systems "Enrico Fermi" 
CXXXIV course (Varenna 1996), edited by F. Mallamace 
and H.E. Stanley (IOS Press, Amsterdam, 1997), p.491. 
M. Nicodemi and A. Coniglio, J. Phys A Lett. 30, L187 
(1996). F. Ricci-Tersenghi, DA. Stariolo and J.J. Aren- 
zon, Phys. Rev. Lett. 84, 4473 (2000). 
[27] A. Coniglio, M. Nicodemi, J. Phys.: Cond. Matt. 12, 
6601 (2000); M. Nicodemi, A. Coniglio, Phys. Rev. Lett. 
82, 916 (1999). 
[28] J.J. Arenzon, J. Phys. A. 32, L107 (1999). 
[29] A. de Candia and A. Coniglio, Phys. Rev. E 65, 16132 



11 



(2002). 

[30] to is measured in Monte Carlo steps (MCS), where 
1 MCS corresponds to N attempts to move a particle 
randomly chosen, and TV = J^. m is the number of par- 
ticles. 

[31] Note that with our approach it is also possible to explore 
low density inherent states in a stationary regime and 
not only the off-equilibrium "glassy regime", as instead 
in Ref. pf. For instance, the frustrated lattice gas model 
at density p — 0.65 is hardly found in an out of equilib- 
rium quasi-stationary state (at any finite value of the 
bath temperature the system quickly reaches the equi- 
librium state). We have considered the model eq.([l]), at 
p = 0.65, and we have performed an usual Monte Carlo 
diffusive dynamics of the model. At a starting time we 
prepared the system in equilibrium with a very high bath 
temperature. Afterwards, it is suddenly cooled at a very 
low Tbath, but the system quickly reaches the equilibrium 
state, and T dyn coincides with T bat h- 



[32] Where not explicitely shown, the error bars of the data 
are more or less equal to the size of symbols in figures. 

[33] J.B. Knight, C.G. Fandrich, C.N. Lau, H.M. Jaeger and 
S.R. Nagel, Phys. Rev. E 51, 3957 (1995). E.R. Nowak, 
J.B. Knight, E. Ben-Naim, H.M. Jaeger and S.R. Nagel, 
Phys. Rev. E 57, 1971 (1998). E.R. Nowak, J.B. Knight, 
M. Povinelli, H.M. Jaeger, and S.R. Nagel, Powder Tech- 
nology 94, 79 (1997). 

[34] In both cases, p = 0.65 and 0.75, the stationary self scat- 
tering functions obtained for finite tap duration To, are 
well fitted by stretched exponentials, eT^ t l T ' . The char- 
acteristic time scale, r, obtained in this way increases 
as Tr decreases, and seems to diverge only in the limit 
T r -> 0. 

[35] In the previous case, p — 0.65, the minimum energy, 
Emi n , is zero. 

[36] A. Coniglio, A. Fierro, and M. Nicodemi (unpublished). 



